High fidelity feed-back assisted parity measurement in circuit QED 



L. Tornberg and G. Johansson 
Department of Microtechnology and Nanoscience - MC2, 
Chalmers University of Technology, SE-41296 Gothenburg, Sweden 
(Dated: March 30, 2010) 

We analyze a two qubit parity measurement based on dispersive read-out in circuit quantum 
electrodynamics. The back-action on the qubits has two qualitatively different contributions. One 
is an unavoidable dephasing in one of the parity subspaces, arising during the transient time of 
switching on the measurement. The other part is a stochastic rotation of the phase in the same 
subspace, which persists during the whole measurement. The latter can be determined from the full 
measurement record, using the method of state estimation. Our main result is that the outcome of 
this phase determination process is independent of the initial state in the state estimation procedure. 
The procedure can thus be used in a measurement situation, where the initial state is unknown. 
We discuss how this feed-back method can be used to achieve a high fidelity parity measurement 
for realistic values of the cavity-qubit coupling strength. Finally, we discuss the robustness of the 
feed-back procedure towards errors in the measurement record. 



I. INTRODUCTION 

In circuit quantum electrodynamics (cQED)JJ^, super- 
conducting qubit J^"^ are coupled to a microwave cavity, 
vifhich allows for high fidelity qubit state controP and 
measurementP^, two ingredients required to achieve a 
scalable quantum information architectural. In the 
dispersive regime, where the cavity and qubits are 
sufficiently detunecf, the state of the qubits shifts the 
resonance frequency of the cavity which can be detected 
by a homodyne measurement of the field state. This 
realizes a joint read-out of the state of the qubits and 
allows for measurement of single, as well as, multi-qubit 
operator^^lSl, Such a resource allows for measurement 
of the state parity which in turn enables deterministic 
entanglement through measurement!^^^^ and imple- 
mentation of error correction piotocol^r^^. A true 
two-qubit parity measurement projects the state of the 
system on either the even parity subspace spanned by 
\gg) and |ee), or on the odd parity subspace spanned 
by \ge) and \eg), simultaneously allowing the observer 
to infer which of the two projections has occurred. The 
measurement is however not allowed to give information 
about the nature of the states within the respective sub- 
spaces as this would destroy the quantum superposition 
of the post-measurement stat^^. 

In the dispersive read-out, the entanglement be- 
tween qubits and field allows the observer to infer the 
parity of the qubit state. As analyzed by Lalumiere et. 
al. in Ref. [13], the same entanglement gives rise to an 
unavoidable back-action causing dephasing within one 
of the parity subspaces. This severely limits the fidelity 
of the measurement and causes the post-measurement 
state to be mixed and thus limited for further use in a 
computational context. 

In this paper, we show that the full measurement 
record from the homodyne current together with quan- 
tum feedback based on state estimatiorPSlSlj can be used 



to significantly reduce the unwanted dephasing caused 
by the measurement. This possibility was previously 
analyzed in a different setup where the homodyne 
detection could be regarded as a strong projective 
measurementPSl. For circuit QED, a weak measurement 
is more realistic, and in the present work we extend the 
analysis to this situation, using a quantum trajectories 
approaclP^. 

II. DISPERSIVE READ-OUT AS A PARITY 
MEASUREMENT 

In this section, we lay out the essential idea behind 
the realization of a two qubit parity measurement us- 
ing dispersive read-out in circuit QED^^. For two 
qubits, the objective of a parity measurement is to dis- 
tinguish between measurement results corresponding to 
states belonging to the two othogonal subspaces 71+ = 
span(|(7(7) , |ee)) and 'H- = spand^e) , |eg)) with corre- 
sponding parities P+ — 1 and P_ = — 1. Furthermore, 
such a measurement is not allowed not distinguish be- 
tween states within the respective subspaces since this 
would destroy superpositions of states in either or 
H- . We consider a cQED system consisting of two super- 
conducting charge qubits coupled to a s ingl e field mode 
of a superconducting stripline resonatop32]. Neglecting 
the infiuence of higher qubit states, the system is well 
described by the Tavis-Cummings HamiltoniarPH 

H = uJra^a + j^'fa^^UY.g^iaa^^+a^af) 

-I- (ae„e " + a' " ), (1) 

where we have set h — 1, ujr is the resonance frequency 
of the cavity mode and e„ is the drive amplitude of the 
measurement signal tuned to frequency w^- The level 
splitting of qubit j is given by ujqj and gj is the coupling 
between the corresponding qubit and the cavity. To re- 
alize a qubit measurement, the system is operated in the 
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dispersive regime where the cavity and qubit frequencies 
are far detuned, A^- = 5j/|Aj| <C 1 where Aj = uir — ujqj. 
In this hmit the system is well described by the second 
order effective HamiltoniarP 

Heff - Ara^a + ^ xA^^a^a + ^^^«±^ai^-) (2) 

I / * I t ^ I T / (1) (2) I (1) (2)\ 

written in the frame where the cavity degrees of freedom 
rotate at the frequency uJm such that — ujr — ujm- The 
residual coupling between cavity and qubit j is described 
by X] = .9|/Aj and J12 = gi52(Ai + A2)/2AiA2 is the 
qubit-qubit coupling mediated by virtual photon^2ll_ 

In the dispersive readout of cQED, the joint state 
of the qubits is inferred from the homodyne signal 
coming from the transmitted microwaves through the 
cavity. Due to the coupling between cavity and qubits, 
the phase and amplitude of the transmitted microwaves 
will depend on the state of the qubits i,j = g,e 
with the field evolving into a coherent state \aij) with 
amplitude a^j which obeys the differential equatiorP^l 

= - «(Ar + Xij)ctij{'t) - (3) 

Here, Xij = (O'lxif^i^^ + X2<^f^ Hi) is the coupling be- 
tween the state \ij) and the cavity. In this case the four 
different cavity states can be interpreted as pointer states 
of the measurement, where the state of the two qubits 
can be inferred from the state of the cavity field. To re- 
alize a parity measurement we must make sure that the 
read-out distinguishes between the two sets of field states 
{aggjttee} and {age,oieg} but not between the ampli- 
tudes within the sets. By choosing the couplings 171 = 52 
and detunings Ai = — A2, the four different amplitudes 
in Eq. ([s]) satisfy 

agg{t) = aee(t), 

Re{age{t)) = -Re{aeg{t)) 

lin{age{t)) = lni{aeg{t)). (4) 

which can be seen in Fig. [T] were the quadrature Q = 
Im(aij(i)) is plotted against the in-phase component of 
the field / — He{aij{t)). In a homodyne measurement 
the field of the cavity is mixed with a local oscillator 
with a relative phase (p defining the measured quadra- 
ture. By choosing (jj = 7r/2, the Q-quadrature is mea- 
sured as indicated in Fig. [T] In this case, the observer 
can distinguish between states of even and odd parity. 
The signal will however not contain information about 
the phase between states {\ge) ,\eg)} and {\gg) , |ee)} as 
required by a parity measurement. Since the measure- 
ment eigenstates are not eigenstates of the qubit-qubit 
coupling (last term in Eq. ([2])) we choose Ai = — A2 to 
make this vanish. 
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FIG. 1: Phase space illustration of the stationary field states 
given in Eq. ([3| with parameters Ar = 0, 6™ = 0.5k, 
Xi ~ — X2 = 1.5k. By choosing the relative phase of the 
local oscillator to be (p = 7r/2, the homodyne read-out mea- 
sures the projection of the field on the Q-quadrature allowing 
the observer to infer the parity of the qubits. 



III. THE MODEL 

A. Effective two qubit stochastic master equation 

The appropriate equation of motion to describe the 
evolution of the system conditioned on continuous homo- 
dyne detection is given by the stochastic master equation 
(SME)24 

dpc ^ CtotPcdt + y/Kr]M[ae^''']pcdW(t), (5) 
where CtotPc is given by 

CtotPc = -i[H^s, Pc] + lipW^-^]Pc + '^V[a^^^]p, 

3 

-f KD[a]pc + nVY Ajcri^'^]pc, (6) 
j 

where 'D[c\p — cpc^ — l/2{c^ cp + pc^c) is a dissipation 
super operator in Lindblad forrrPni with the pure relax- 
ation and dephasing rates of qubit j given by 71-,- and 
70j respectively. The cavity damping rate is given by n 
and the last term of Eq. ([6| describes the effect of Purcell 
dampinJSH. The measurement back-action on the system 
is described by the measurement super-operator 

M[c]pc ^ cpc + PcC" - {c + C^) Pc, (7) 
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where (•) = tr(-pc) and dW{t) is defined as a Wiener 
increment completely characterized by its mean and 
varianc^^ni 



^dW{t)] = 0, 
E[dW{t)^] = dt. 



(8) 



Here, E[-] denotes the ensemble average taken over differ- 
ent realizations of the noise process W{t). The homodyne 
current is given by 



j{t)dt = ^{ae-'f + a}e'^)dt + dW{t), 



E 



(10) 



dpc 
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+ E X*:'>*j (I™("ii"feO + iR-e(aijafc/))ni-,Pcnfe;(ii 

ij,kl 

+ y/J^M[Il„e-"t']p^dW{t), 

where XijM = Xij~Xki- The fifth term represents the de- 
phasing induced by the measurement and the sixth term 
gives the AC Stark shift. The measurement operator Ha 
is given by 



By choosing system parameters appropriately, the homo- 
dyne measurement along with single qubit rotations can 
thus be tuned to measure single, as well as two qubit op- 
erators as discussed in Refs. [7 ^9113] . The current can 
be expressed in terms of qubit operators as 

j{t)dt = ^]{Uae-''l' + Ule">')dt + dW{t), (14) 

which in the case of (p = tt/2 reduces to 



(9) 

where rj is the efficiency at which the photons are 
detected. 

The aim of this section is to trace out the cavity 
degrees of freedom from Eq. ^ to obtain an effective 
master equation for the two qubits only. As shown in 
Refs. |13|27j this can be done in the limit jij <^ k, which 
is easily satisfied with current Purcell limited qubit^^H, 
In a frame of reference given by the transformation 



jij{t)dt = 2y/KT]lT[i{aij)dt + dW{t), 



the photon population in the cavity is essentially zero, 
which allows one to trace out the cavity degrees of free- 
dom. Here, D[a] — exp(Q;a^ — a* a) is the displacement 
operator of the fielcP® and H^ = \ij) {ij\ are projection 
operators onto the respective basis state of the two qubit 
Hilbert space. This gives an effective master equation for 
the qubit degrees of freedom 



(15) 



for the basis states With our choice of parameters, 
Ini(Q;gg) — Im(Q!ee) and Im(age) = Im(Q;eg) such that the 
measurement does not distinguish between states within 
and "H- as required for a parity measurement. 



IV. MEASUREMENT INDUCED DEPHASING 

With the field amplitudes given in Eq. Q and the LO- 
phase given by </) = 7r/2 the measurement operator can 
be rewritten as 



(16) 



where 



n+ = Igg) {ggl + \ee) (ee| , 
n- = |.ge) {ge\ + \eg} {eg\ , 



(17) 



are projection operators onto and and /3 is the 
difference field 



f3 = Im(agg) - Im(Q;ge) 



(18) 



pdff he first term of Eq. ( 16 ) is associated with information 



gain and simply expresses the fact that the parity mea- 
surement tends to localize an initial state in one of the 
subspaces or "H- at a rate given by 



(11) 



(19) 



H, 



(12) 



Izicri^^ + liz^f^ + 711 



which can be taken as a definition of the measurement 
rate. The second term represents a stochastic phase 
between the states \ge) and \eg) accumulated during the 
measurement. This, together with the dephasing term 
oc Im(a*jafc;) given in Eq. (11), represents the unwanted 
back-action of the measurement which we now analyze. 

Because of the non-linear nature of conditional state 



evolution, it is not possible to analytically solve Eq. ( 11 ) 



with the amplitudes 

^zz 
Izl 
llz 

711 



^99 ^ge ^ ^eg i '-^ee? 
~^99 ^ge ^ea "~t~ Q^g 



in the general case. However, in the case when the state 
has been projected onto 'H± the non-hnear part vanishes 
and Eq. (11) reduces to 



zgg I i^ge ^eg ^ee-) 



dpc{t) 



(13) 



^-V[Y.z\pc(t)dt - iy^^lCpc{t)dW{t), 

(20) 
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where E2 = Ilge - Ileg, JCpc = P^,^], Tg^^eg = 

4xlm(a*gaeg) and = Re^(age)- Since the focus of the 
analysis is on the back action of the measurement we have 
neglected the effect of pure relaxation and dephasing as 
well as Purcell damping. Furthermore, we disregard the 
deterministic rotation of the phase oc Re(a*jQ!fei) given in 
Eq. (11) since this can be u ndone regardless of the mea- 
surement outcome. Eq. (20) can be solved analytically^^, 
yielding 



Pc{t) 



exp 



rge,eg(s)X>[S2] + Krj^l{s)}C'^ ds 



-2.yo 

i / ^/Kr]n{s)dW[s)K\p[0) 



exp 



{Vge,eg{s) - 2K.riVl{s))'D[Y,^]ds 



i / ^/Krin{s)dW{s)Jc\p{0), 



(21) 



where we have used the fact that IC^p = —2T>[T,z]p- 
Given this, we see two effects of the measurement. First, 



the fourth line in Eq. (21 ) gives the dephasing rate Tc for 
a single trajectory 



(22) 



The second effect is given by the last term in Eq. (21) 



which gives rise to a stochastic phase between \ge) and 
\eg). Since this phase varies between different measure- 
ments the ensemble averaged state is mixed with the cor- 
responding dephasing rate given by Fge.eg- As a measure 
of the read-out fidelity we take the ratio between the 
measurement and dephasing rate as t — ?> 00 



(23) 



which shows that the dispersive read-out only works as 
a perfect parity measurement in the large coupling limit 
X S> K. This limit may however be experimentally hard 
to reachl'^ and it is therefore desirable to seek alternative 
methods to improve upon this ratio. Considering the 
same ratio but for a single trajectory gives the result 



L |t— >0 



(1 - 77)k2 



(24) 



which vanishes in the case of a perfect measurement 
ry = 1. This can be understood by looking at the two 
terms in Eq. (22). In Fig. [2j we see that for a perfect 
measurement, Fc — > as t — !■ cx), which is due to the 
fact that the measurement current contains information 
about the photon fluctuations. As i — >■ 00, the observer 
can track the photon fluctuations by looking at the ho- 
modyne current such that no information is lost into the 
environment. Given this information, the observer can in 
principle undo the stochastic phase with a resulting pure 
ensemble averaged state. From Fig. [2j we also see that. 



for the dispersive measurement, there is initial dephasing 
in the system not caught when taking the limit t ^ 00. 
This is given by 



lim 



Tc{s)ds 



(^2 -f 16x2)2 ^ 



7? = 1, (25) 



which simply reflects the fact that the homodyne 
measurement is not able to record information about 
the initial photon fluctuations due to the finite response 
time of the cavity. Because of this, inevitable dephasing 
occurs with the final state being mixed even for a single 
trajectory. The dephasing in Eq. (25) is a decreasing 



function of x since the number of photons in the cavity 
interacting with the qubits decreases as the coupling to 
the qubits shifts the cavity frequency. Moreover, the 
dephasing decreases with increasing k. This is simply 
due to the fact that a faster response time of the cavity 
allows for a more efficient measurement of the initial 
photon fluctuations. 

From Eq. (24 1, we note that the measurement fi- 
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FIG. 2: (Color online) The dephasing rates Fge.es (solid blue 
line) and 2K.riQ (dashed red line). Initially, the system is prone 
to dephasing wth Fge.eg > 2Hirj^l. This reflects the fact that 
initial fluctuations in the photon number can not be captured 



by the homodyne measurement. As t — >■ 00, Pg 



with all the unwanted back action given by a stochastic phase 
between the basis states \ge} and \eg}. The parameters are 
= K, xi = -X2 = 3k and 77 = 1. 

delity of a single trajectory is far better than that given 
for the ensemble averaged state in Eq. (23). This is 



because the stochastic phase does not contribute to the 
dephasing for a single trajectory. From this it is clear 
that the purity of the ensemble averaged state can be 
drastically increased if the observer is able to undo the 
stochastic phase after each measurement shot. In this 
case the dispersive read-out approaches a true parity 
measurement even in the case of moderate coupling 
between qubits and cavity. In the next section we 
present such a scheme based on state estimation and 
feed-back. We show that this can be done without 
knowledge of the initial state. 
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V. INCREASING MEASUREMENT FIDELITY 
BY STATE ESTIMATION 

Information about the stochastic phase discussed in 
Sec. llVl can not be extracted from the measurement cur- 
rent directly. We may however use the recorded current 
to extract the noise process W{t) for a given trajectory. 
This can be used as an inpu t to Eg. (11) to calculate 
an estimate Pc(i) of the stat d^^ l ^^ . Unfortunately, this 
procedure requires knowledge of the initial state Pc{0) 
which is typically unknown in a measurement situation. 
We now show that no information about the initial state 
is needed if we only are interested in the accumulated 
stochastic phase. 

Given an unknown initial state, the observer will 
make a wrong estimate dW{t) of the true noise process 

dW{t) = dW{t) (26) 
+ V^tr [(n^e-''^ + nt„e'*) (p,(t) - p,(i))] dt, 

where Pc{t) is the true state of the system and Pc{t) the 
guessed estimate. We write the true state as a mix of 
odd and even parity states 



state is a mix or superposition between even and odd 



(27) 



where p+{t) G 7i+ and p-{t) G The probability 

p{t) G {0, 1} describes the conditional evolution towards 
one of the respective subspaces. We choose the estimated 
state to initially reside in the 7i- subspace such that the 
subsequent state can be written 



p,{t)^p_{t). 



(28) 



Considering the measurement part of Eq. (11) only, we 
can write three coupled equations for the \ge) ( eg\ matrix 
element of pc, pc and p{t) 

dpge,eg = -2y/KJ^[l3p + iRe{age)) Pge,egdW {t) , 
dpge.eg = -i2y^{2y/Krjl3pdt + dW {t))Re{age)pge,eg, 

dp{t) = -2^Pp{l~p)dW{t). (29) 



Given the solution for p{t), we can calculate the phases 

4> = arg(/cige,eg) and = arg(/5ge,eg) of the true and esti- 
mated state respectively 

0(i) = ^{t) = -Anri [ /3{s)Re{ageis))p{s)dt 
Jo 

- 2^ I Re{age{s))dW{t). (30) 
Jo 

Hence, we see that the accumulated phase of the esti- 
mated state equals that of the true state! This is the 
main result of this paper. This is consistent with the 
result of Ref. [20], where a strong projective homodyne 
measurement is considered. This can be seen as the lim- 
iting case when k — > oo of our analysis. Since the con- 
ditional evolution does not depend on whether the true 



parity states, the main result in Eq. (30) holds equally 



well when the true state is an arbitrary superposition 
state. Hence, the method of undoing the accumulated 
phase can be applied not only in the context of e.g. en- 
tanglement generation as considered in Ref. [13] but also 
in a real measurement situation where the initial state 
is unknown. This will considerably reduce the unwanted 
dephasing casued by the measurement back action. 



VI. RESULTS 

To quantify how much the method of state estimation 
and subsequent subtraction of the accumulated phase in- 
creases the purity of the post-measurement state we per- 



form numerical simulations of Eq. (111. This is necessary 



since the full dynamics of Eq. ( 11 ) can not be solved an- 



alytically. To distinguish the odd and even parity states 
we use the integrated current 



s{t) 



jis)ds, 



(31) 



and assign the post measurement states to 'H+ or if 
s{t) satisfies s{t) > Stu or s{i) < Stu where Stu is a thresh- 
old given by the average current Sth = Jq (Ini(a3g) + 
hii(age))ds. To quantify the fidelity of the parity mea- 
surement and the estimation scheme we define the aver- 
age concurrence C = C(E[pc]), purity P = tr((E[pc])^) 
and fidelity 

-^ n+(ij+\ E+ [p j + n_ ( I E_ [p,] \^_) ^g^) 

where n^/ri- is the number of trajectories assigned to 
measurement outcome +/— and E-|-[pc] is the ensemble 
average over the respective states. The states \iIj±) are 
even and odd parity states with the initial state given by 



IV') =cos(6i)|V+)+sin(6i)|V-) 



(33) 



The initial estimated state is correspondingly given by 



ip- ) where 



need not equal We 

plot C, P and F in Fig.|3l Fig. H] and Fig. [s] for different 
values oi 9. In each figure the inset shows the value 
of the respective quantity at the final time with and 
without the subtraction of the accumulated phase. It is 
clear that the state estimation allows us to subtract the 
stochastic phase in each measurement run, resulting in a 
drastic improvement of the state fidelity across all three 
measures. 



A. Imperfect measurement record 

We note that our estimation scheme relies on the as- 
sumption that we can make a perfect record of the homo- 
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FIG. 3: (Color online) Average concurrence C as a function 
of time plotted for five initial states on the form given in 
Eq. (27 1. The angle 6 is in the range 9 G {0, 7r/4} (increasing 
from top to bottom). The inset shows C at the final time 
of the measurement, with (blue o) and without (red o) the 
phase subtracted. The parameters are given by e™, = k, xi = 
-X2 = 3k, 31 = -32 = 100k, 7ij = = and r; = 1. 
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FIG. 5: (Color online) Average fidelity as a function of 
time plotted for the same initial states as in Fig. |3] The inset 
shows F at the final time of the measurement, with (blue o) 
and without (red o) the phase subtracted. The parameters 
are the same as in Fig. |3] 
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FIG. 4: (Color online) Average purity P as a function of time 
plotted for the same initial states as in Fig. |3] The inset 
shows P at the final time of the measurement, with (blue o) 
and without (red o) the phase subtracted. The parameters 
are the same as in Fig. [3] 



dyne current, which is impossible in the case of a detector 
with limited bandwidth. We therefore consider the situ- 
ation where the estimated state is based on the value of 
the low pass filtered current 



iip(i) 



j{s)ds, 



(34) 



where r gives a measure of the deviation of Jip(t) 
from j{t). In this case the value of the estimated 
phase will differ from the true value and the averaging 
over measurement results will reduce the purity of 
the state. In Fig. |6] we plot the imaginary and real 
part of Pge,eg of the post-measurement density ma- 
trix for each trajectory. The initial state is given by 
IV') = l/v/2(|ge) -hexp(i7r/4) \eg)). We see that the state 
estimation scheme successfully recovers the true phase 
of the initial state even in the case where we do not have 
access to a perfect measurement record. As r increases 
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FIG. 6: (Color online) In a)-e), we plot Im(pge,eg) of the post 
measurement states against Re(pge,es). The initial state is 
given by \ip) = l/y/2{\ge) + exp(i7r/4) \ eg)). As r increases 
the estimation of the accumulated phase becomes worse and 
the phase of the state is distributed around — ~tt/4. In 
f), the average purity is plotted as function of r. Even for 
T = 50dt, the purity does not reach the corresponding value 
obtained in the case without state estimation and feedback. 
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the fidelity goes down as expected, but does not reach 
the value obtained without state estimation even for 
T — 50dt. We can understand this robustness from the 
observation that the values of Pge,eg are distributed on 
a circle arc which make the decrease of purity quadratic 
in the phase uncertainty. This can also be seen in the 
initial slope in Fig. [6j 



VII. CONCLUSION 

In conclusion, we have performed an analysis of a two 
qubit parity measurement based on dispersive read-out 
in cQED. In particular, we analyzed the the back-action 
on the two qubits and found two qualitatively different 
contributions. One is an unavoidable dephasing in one 
of the parity subspaces, arising during the transient time 
of switching on the measurement. This dephasing occurs 
in the process of entangling the state of the driven cav- 
ity with the two qubits. The other part is a stochastic 



rotation of the phase in the same subspace, which per- 
sists during the whole measurement. We discussed how 
the latter can be determined from the full measurement 
time trace, using the method of state estimation. Quite 
surprisingly, we found that the outcome of this phase de- 
termination process is independent of the initial state in 
the state estimation procedure, making it useful in the 
situation of a parity measurement, where the initial state 
by definition is unknown. Finally, we show that the feed- 
back process is rather robust towards imperfections in 
the measurement record. This analysis opens up for re- 
alizing high-fidelity parity measurements in circuit QED, 
using realistic values of the coupling strength between the 
qubits and the cavity. 
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